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ABSTRACT 

The Hertz and SCUBA polarimeters, working at 350 yum and 850 /im respec- 
tively, have measured the polarized emission in scores of Galactic clouds. Of the 
clouds in each dataset, 17 were mapped by both instruments with good polariza- 
tion signal-to-noise ratios. We present maps of each of these 17 clouds comparing 
the dual-wavelength polarization amplitudes and position angles at the same spa- 
tial locations. In total number of clouds compared, this is a four-fold increase 
over previous work. Across the entire data-set real position angle differences are 
seen between wavelengths. While the distribution of 0(850) — 0(350) is centered 
near zero (near-equal angles), 64% of data points with high polarization signal- 
to-noise (P > 3ap) have |0(85O) - 0(350) | > 10°. Of those data with small 
changes in position angle (< 10°) the median ratio of the polarization ampli- 
tudes is P(850)/P(350) = 1.7 ± 0.6. This value is consistent with previous work 
performed on smaller samples and models which require mixtures of different 
grain properties and polarization efficiencies. Along with the polarization data 
we have also compiled the intensity data at both wavelengths; we find a trend 
of decreasing polarization with increasing 850-to-350 /im intensity ratio. All the 
polarization and intensity data presented here (1699 points in total) are available 
in electronic format. 



Subject headings: dust, extinction — ISM: clouds — polarization — submil- 
limeter: ISM 
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Introduction 



Observations of polarized radiation in the interstellar medium at optical through mil- 
limeter wavelengt hs have been attributed to extiriction by, and emission from, interstellar 
dust grains (e.g., iHiltnerl Il949l Il95ll : iHildebrandl Il988[ ). In order to generate a net po- 



larization the grains must be aspherical and exhibit a relative net alignment of their axes 
with one another and with the interstellar magnetic - field, typically the sh ortest grain axis 
is parallel to the field (e.g., iDavis fc Greensteinlll95ll : lLazarianll2003l . 120071 ). At near-optical 
wavelengths this polarizing dust screen causes the dichroic extinction of background starlight 
with respect to the grain axes' different cross-sections. At far-infrared and longer wavelengths 
dominated by grain emission rather than extinction the polarization results from the axes' 
different emission cross-sections. Due to the necessary role of magnetic fields in aligning 
the grains, polarization observations have been us e d primarily to s t udy interstellar magnetic 



fields fe.g.. iHildebrand 



20071 : IPereyra fc Magalhaes 



19881: iFosalba et al.l l2002l : ICrutcherl 12004 : ICurran fc Chrvsostomou 
2007h . cifically, the magnetic field morphology (projected 



onto the plane-of-the-sky) is inferred from the polarization position angles. (The field is 
typically parallel to the polarization angle in the case of extinction and perpendicular in 
the case of emission.) However, the physical properties of the grains themselves and their 
interaction with the field are no less important than the field itself; many of these proper- 
ties can be inferred from the polarization amplitude. For example, spectropolarimetry of 
background- starlight at near-optical wave lengths has been used to measure the shapes of 
dust grains (IHildebrand &: Dragovanlll995), ma ke tests of grain alignment mechanisms (e.g.. 



Whittet et al. 



2008 



Kim fc Martin 



Andersson &: Potterl 120071 ) . place limits on the size of aligned grains 



(e.g., 

spectral lines (e.g.. IWhittetl 120041 ). 



19951). and m easure the composition of the aligned grains via polarized 



Polarimetry at submillimeter wavelengths was initially driven by the desire to study the 
magnetic-field morphology of interstellar clouds. As such, the amplitude of the polarization 
(typically ~1-10%), and any wavelength dependence, was mostly secondary to measurements 
of the polarization position angle. Most studies of the polarization spectrum in the far- 
infrared and submillimeter have thus relied on observations where the choice of wavelength 
was made according to the availability of atmospheric observing windows, not with the 
specific goal of studying any spectral variation itself. This mode of operation h as resulted in 



a num ber of da tabases at a few spe cific wavelengt hs including 60 and 10 /xm ( iDotson et al. 



2000h . 350 /xm JPotson et allboioh. and 85 /im jMatthews et al.l 120091 ). Using a subset of 
the available data, IHildebrand et al.l ( 119991 ) showed that the polarization spec trum across 
these four wavelengths had a minimum near 350 /im (see also IHildebrandl 1200 ll ). 



The spectral structure observed in near-optical continuum polarimetry (i.e., the "Serkowski 



- 3 - 



law" ; ISerkowski. Mathewson. fc FordI Il975[ ) is the result of a combination of properties: a) 



interstellar dust grains have typical radii a ~ 0.1-1 yum (e.g., iMathis. Rumpl. &: Nordsieck 
19771 ): b) the larger grains are more efficiently aligned than the smaller grains; and c) 
the grain sizes prob ed are of the same order as the observing wavelengths A ~ a (e.g., 
Kim fc Martini 119951 ). On the other hand, at wavelengths in and beyond the far- infrared 
(A > 50 /zm) all the above properties are the same save for the fact that the grains are com- 
paratively small (A ^ a). In this case one expects no variation in t he polarization spectruni . 
Therefore, any model explaining the spectral structure observed by iHildebrand et al.l (119991 ) 
requires multiple dust grain populations in which there is a correlation between the efficiency 
with which the grains are aligned and other properties related to the emitted radiation (i.e., 
grain size, temperature, emissivity). We describe some simple models in Section |H 

The initial studies of Hildebrand et al. Jl999h did not have sufficient data to test such 
models. Vaillancourt ( 20021 ) and Vaillancourt et al. ( 2008 ) extended these datasets shghtly, 
performing clo ud-by-cloud comp arisons as well as point-by-point spectral comparisons within 
clouds. While IHildebrand et al.r s original result held-up under this more detailed analysis, 
the later work was still limited to a small number of molecular clouds. The recent compilation 
of large re-reduced datasets at 350 /im and 850 /xm presents the opportunity to further extend 
the sample to point-by-point comparisons in additional Galactic molecular clouds. In this 
work we present a comparison of submillimeter polarization data at these wavelengths in a 
total of 17 clouds. 

We present polarization maps of these 17 objects comparing the polarization amplitude 
and position angle at the two wavelengths. Section [3] highlights differences between the two 
wavelengths in both angle and amplitude. Changes in the angle may help disentangle the 
magnetic field morphology alon g the line of sight or e xtend maps to regions not observ- 



able at some wave lengths (e.g., ISchleuning et al.ll2000l : iKandori et al.ll2007l : iLi et al.l 12009 



VaillancourtI |2012| ) . but we do not elaborate on the position angle data presented here. In 
Section [3] we also compare the 850-to-350 yum polarization ratio on a point-by-point basis. 
In Section HI we compare both the absolute polarization magnitudes, as well as the 850-to- 
350 /xm polarization ratio, with the 850-to-350 /xm intensity ratio, compare it to previous 
work, and briefly discuss grain alignment models consistent with the data. All the polariza- 
tion and intensity data presented in this paper are available in machine readable tables in 
the electronic version of the journal. 
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2. Data Processing 



From ~ 1995 to 2005 independent campaigns to map the polarization at 350 /xm and 
850 fim were carried out by the Hertz polarimeter at the Caltech Submilhmeter Observatory 
(CSO) and the SCUBA polarimeter at the James Clerk Maxwell Telescope (JCMT), both 
on Mauna Kea, Hawaii. The Hertz passband was chosen to match the 350 fim atmospheric 
window while SCUBA's polarimeter operated primarily at 850 /im (Figured]). 



2.1. Spatial Resolution and Map Sampling 



The H ertz instrurnent, i t s observing strategy, and data analysis are described in deta il 
elsewhere Jxirbv et al.l bood : lOowell et all Il998l : ISchleuning et al.l Il997l : IPlatt et all Il99lh . 
Here we briefly review the aspects relevant to the present work. Hertz incorporates two 
separate detector arrays which simultaneously measure the two orthogonal modes of linear 
polarization, modulated by a cryogenic half- wave plate (HWP). The 6x6 pixel^ arrays have 
pixel center-to-center spacings of 17'.'8 and a beamsize of approximately 20" full-width at 
half-maximum (FWHM). The observing strategy involves rotating the instrument to follow 
the sky- rotation throughout the night, as well as steps of order the array-size to map areas 
larger than the 2' x 2' field-of-view. The rotation allowed a single pixel to continue observing 
the same patch of sky throughout the night. Additionally, step-sizes were typically chosen 
to be an integer-number of pixels; only rarely were maps made with samples spaced more 
closely than the 17'.'8 pixel pitch. As a result, these data do not meet the Nyquist criterion for 
a fully-sampled map and we have, therefore, made no attempt to gener ate polarizatio n maps 
at increased spatial sampling. All the polarization data reported by iDotson et al.l (120101 ) 
and used here maintain a spat ial resolutiori of 20" . (Measured beam profiles are given here 
in Figure [2] and in Figure 3 of iDowell et al.l ( 119981 ).) 



The SCUBA camera and polarimeter is described in detail elsewhere ( iGreaves et al. 



and polarime 


ter is 


Holland et al. 


1999) 



consists of a single detector array with 37 pixels arranged on a hexagonal grid with a 2'.3 
field-of-view; the polarization is measured by inserting a warm wire-grid, modulated by a 
stepped HWP. Fully-sampled polarimetric and photometric maps are generated by "jiggle- 
mapping" which moves the array by sub-pixel steps. The data analysis involves combining 
the individual jiggle maps and re-sampling them onto an output grid with a 6" pitch. Gen- 
eration of fully-sampled maps in this manner also alleviates the need to follow sky-rotation 
throughout the night. While the intrinsic SCUBA beam size at 850 fim is close to the 14" 
diffraction limit of the JCMT, this map-making process results in an effective beam-size 
closer to 20" (Figure [2]). 
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Fig. 1. — The Hertz 350 /xm and SCUBA 850 /iin passbands (solid curves) are shown along 
with the typical atmospheric transmission on Mauna Kea (dashed line). The atmospheric 
transmission is calculated for 1mm of precipitable water vapor using the CSO's web-based 
interface (http: //www.submm.caltech.edu/cso/weather/atplot.shtml ). 
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Fig. 2. — The Hertz 350 fim (circles and dashed hne) and SCUBA-pol 850 /xm (diamonds 
and dash-dotted line) beams were measured on Uranus in September 1997 and July 1998, 
respectively. In both cases the disk's angular diameter was 3'.'7. These data represent the 
effective beam-widths after all data processing steps. Gaussian beams with full-widths at 
half-maximum ranging from 17" - 21" are shown for comparison (gray). Error-bars on the 
Hertz data indicate the variation of the signal within an annulus centered at each radius 
(IDowell et al.l 119981 ). Error-bars on the SCUBA-pol points were calculated from intensity 
uncertainties at each data pixel and are smaller than the plotted data points. 



- 7- 



As noted above, the resultant spatial resolution for both Hertz at 350 fim and SCUBA- 
pol at 850 /im is ~ 20". We refer to this resolution as the effective beam- width, which 
should not be confused with the diffraction limited resolution of the respective telescopes. 
Direct comparisons of the two beams measured on Uranus are shown in Figure [H (The size 
of Uranus at the time of each measurement was 3'.'7.) Gaussian- fits to these profiles yield 
19'.'5±0'.'1 for SCUBA-pol and 19" ±2" for Hertz; the reported statistical uncertainties follow 
from a formal non-linear fi t to the gauss i an pro file. These beamsizes are consistent with the 



measu rements reported by iDowell et al.l (119981 ) for Hertz (20" ± 2") and iDi Francesco et al. 



( 120081 ) for SCUBA (a primary beam of 19"5). Given this beam similarity we have made no 



correction for different spatial resolutions when comparing these data sets. 

The Hertz data are undersampled. However, since the SCUBA data are fully-sampled, 
there is sufficient information to estimate the SCUBA intensity and polarization at the same 
sky locations of the Hertz data. This is accomplis h ed by reducing the SCUBA-pol data in 



the same manner as presented in [Matthews et al.l (120091 ) but choosing to output the data 



to grids and map-center locations which match the Hertz dataset. Table [T] lists the objects 
observed by both Hertz and SCUBA-pol at 350 and 850 /im; Table 2 (in the electronic 
version only) gives a more complete list of the locations within each of the clouds. Table 2 
also includes data for all points at both wavelengths for the polarization magnitudes and their 
ratio P(850)/P(350), position angles and their difference (/)(850) — 0(350), intensity values 
and their ratio F(850)/F(350), and uncertainties on all values. All polarization magnitudes 
in Table 2 have been corrected for positive bias (Section 12. 2p . The best estimates of those 
values is sometimes zero; as a result the ratio P(850)/P(350) is reported as Nan for cases in 
which P(850)=P(350)=0, Inf for cases in which only P(350) = 0, and equal to zero when 
only P(850) = 0. 



Table 1. Object Summary and Polarization Ratios 



data satisfying Pc > data satisfying Pc > 3ap ...also satisfying lA^I'^ < 10° 



Source 


TotaP 


Number^ 


Median 


MAD 


Number*^ 


Median 


MAD 


Number^ 


Median 


MAD 


W3 


91 


58 


2.2 


1.1 


14 


2.5 


0.8 


4 


2.0 


0.3 


NGC 1333 


154 


70 


2.1 


1.4 


3 


3.5 


1.8 


1 


6.5 




OMC-1 


240 


203 


1.6 


0.7 


136 


1.6 


0.6 


65 


1.6 


0.5 


OMC-2 


90 


54 


1.9 


0.9 


10 


2.4 


0.6 


2 


2.6 


0.3 


OMC-3 


99 


82 


2.3 


1.1 


35 


2.0 


0.8 


15 


1.6 


0.6 


NGC 2024 


104 


72 


2.3 


0.9 


19 


2.8 


0.6 


5 


2.2 


0.2 


NGC 2068 LBS 10 


62 


45 


1.7 


0.6 


23 


1.9 


0.5 


7 


1.9 


0.8 


NGC zObo LhSD 17 


bz 


z4 


l.D 


O.D 
















NGC 2071 


50 


18 


4.0 


2.0 


1 


2.4 











Men R2 


76 


58 


3.2 


1.4 


22 


2.8 


0.9 


1 


1.0 




MonOB112'i 


66 


47 


3.7 


2.1 


11 


4.8 


1.6 


1 


9.8 




pOph 


100 


76 


1.7 


0.8 


25 


1.3 


0.4 


7 


1.3 


0.5 


IRAS 16293-2422 


63 


18 


5.2 


1.9 
















NGC 6334A 


59 


37 


2.7 


1.6 


5 


3.0 


0.5 


2 


4.3 


1.8 


W49A 


55 


46 


2.6 


1.4 


23 


2.4 


0.9 


7 


2.7 


0.8 


W51 A (G49.5-0.4) 


112 


74 


3.6 


2.2 


15 


4.5 


2.6 


2 


6.0 


4.2 


DR2P 


216 


142 


1.6 


0.6 


56 


1.6 


0.4 


22 


1.6 


0.4 


DR21 (Main) 


100 


72 


1.7 


0.6 


36 


1.6 


0.5 


14 


1.6 


0.4 


All 


1699 


1124 


2.1 


1.0 


398 


1.9 


0.7 


141 


1.7 


0.6 



Note. — The median polarization ratios (P[850]/P[350]), the median absolute deviation of their distribution (MAD; eq. [2]), 
and the number of data points in each sample are shown here. The column labeled "Pc > 0" indicates regions where the measured 
data satisfy Pm/fp > \/2 at both wavelengths (see Section [2^2]) . Columns labeled "Pc > 3(Tp" indicate only data points satisfying 
that signal-to-noise criterion at both wavelengths; these values are calculated after applying the de-biasing technique discussed in 
Section [2.2[ The columns labeled "also \A(f)\ < 10°" satisfy both th e latter P/cp constraint as well as the additional A<j) criterion.'^ 



Source coordinates can be found in Table 1 of lDotson et al.l (|201Ci ). 



'^Total Number of points common to both the 350 and 850 /im data sets. 

^Number of data points from the "Total" column satisfying the criteria above. 

'^Acj) = 0(850) — 0(350), where 4> refers to the polarization position angle at the noted wavelength. 

<^IRAS 06382-F0939 

^All data in DR21, including DR21 (Main). 



CO 
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2.2. Positive Polarization Bias 



By definition, the polarization amplitude is a positive-definite quantity. As a result, a 
noisy measurement of a truly zero-polarization source will result in a mean positive polar- 
ization. While there is no exact analytical method to correct for this bias many authors use 
the formula: Pc ~ {P^ — cr^)^/^, where Pc is the bias-corrected pola rization, Pm i s the mea- 
sured polarization, and a„ is the measu red polarization uncertainty. IVaillancourtI (l2006t also 
Simmons fc StewartI Il985l : iQuinnI |2012| ) showed this wa s a good estimat o r whe n P^ > 3crp 
but that Pc = was a better estimator if Pm/cTp < a/2. iMatthews et al.l (120091) ap phed the 
formula for high signal-to-noise to all their 850 fim data while iDotson et al.l ( 120101 ) applied 
no corrections to their 350 /im data. For the data presented in this work we set P^. = in 
cases where Pm/o'p < a/2 and apply the above formula for data with Pm/c^p > V^- When 
computing signal-to-noise cuts on P/ cxp in Table [1] and the following sections we use the 
corrected values Pc as described above. The conclusions drawn in Sections [3] and H] use only 
the Pc > 3crp sample and are therefore not effected by the fact that the correction does not 
provide the best estimate within the regime y/2 < P/Op < 3. 

The corrections on P have no effect on estimates of the polarization position angle; 
that is in the sense that there is no bias in the angle estimate as there is for the polar- 
ization amplitude. However, there is necessarily an effect on the angle uncertainty (e.g., 
Naghizadeh-Khouei &: Clarkelll993l ). For example, in cases where the Pm/cp < V2 the best 
estimate is Pc = and thus any angle measurement is meaningless. We have made no at- 
tempt to estimate "correc ted" angle uncertainties in Table 2. Such corrections are small for 
high signal-to-noise data (INaghizadeh-Khouei fc Clarkd 119931 ) and therefore they will have 
little effect on our analysis and discussions below which use only data with Pc > 3crp. 



3. Comparison 

Maps comparing the 350 /xm and 850 /xm polarization data in each region are shown in 
the Appendix. Here we concentrate on quantitative comparisons of the polarization magni- 
tudes and angles between the two different wavelengths. 

The uncertainties on individual polarization measurements and on combined quantities 
like their angle differences and ratios can be quite large. In the measurements below we 
concern ourselves with the question of whether the observed distributions arise solely from 
measurement uncertainties or also have significant contributions from intrinsic variations 
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across the clouds and/or locations within those clouds. To do this we define the reduced-x 

N 



2. 



i=l * 

where Xm is the median value of the samples Xi, cTi is the measurement uncertainty on the 
quantity Xj, and N is the total number of data points in the sample. 



3.1. Polarization Angles 

One of the most obvious aspects of the polarizations presented in the maps is the 
agreement, or disagreement, of the position angles between the two wavelengths. Figure [3] 
compares these angles across the entire data-set and within three clouds with the largest 
number of data points. All the distributions peak near angle differences of zero degrees, that 
is A0 = 0(850) — 0(350) = 0°. For all non-zero polarization data (i.e., data where both 
P(350) > and P(850) > 0) the median angle difference is 1° with a standard deviation 
of 39°. This agreement is stronger (i.e., has a smaller deviation) if we limit ourselves to 
only points with P > 3ap (398 points in the entire data set). In this case the median angle 
difference is 4° with a standard deviation of 28°. 

To rule-out the possibility that the width of the angle distribution is strongly dependent 
on the measurement uncertainties we calculate the Xr value of the angle differences following 
equation ([1]). If most data points were consistent with the median angle difference within 
their uncertainties then we would expect Xr ~ 1 (especially for A — 1 = 397 degrees-of- 
freedom for the entire 3a data-set). However, for the complete data set with P > 3ap we 
find Xr = 14. From this we conclude the distribution's width is intrinsic and not a result 
solely of data uncertainties. 

Most of the clouds in our sample have insufficient data to perform this analysis separately 
on individual clouds. Exceptions to this point are OMC-1, OMC-3, and DR21, whose angle 
distributions are also shown in Figure [31 For data satisfying the P > 3ap criterion in those 
three clouds the median angle differences and standard deviations are 3° ± 25°, 4° ± 19°, and 
10° ± 22°, with xl values of 18, 5, and 7, respectively. Therefore, as was observed for the 
entire data set above, the width of the angle-difference distribution (A0) in these individual 
clouds is real, in the sense that they are not a result solely of the measurement uncertainties. 

Lastly, we should note that the position angle rotations observed in these clouds are 
unlikely to be the result of Faraday rotation. In typical interstellar cloud conditions, at 
these wavelengths, Faraday rotation is generally much srn aller than that observed here (e.g.. 



Schleuning et all l2000l : Matthews. Wilson, fc Fiegell200lh 
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Fig. 3. — Histograms of the 850-to-350 /iin polarization angle differences, (psbo ~ 0350; for o,) 
all clouds in this work, b) OMC-1, c) OMC-3, and d) DR21. In each panel histograms are 
shown for all non-zero polarization data (open histograms) and only data satisfying P > 3crp 
(gray histograms). The bin- widths for (a) are 5°, all others are 10°. The total number of 
data points in each sample is given in Table [H 
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3.2. Polarization Ratio 

Table [1] lists the total number of locations where measurements were made at both 
wavelengths. For the best comparisons we typically choose to include only data satisfying 
the signal-to-noise criterion P > 3ap] this criterion is applied after the de-biasing correction 
discussed in Section 12.21 Figure |l]a shows the distribution of these points. This figure shows 
all data satisfying the Sap criterion, including outliers as high as P(850)/P(350) = 31; the 
inset concentrates on the main peak. 

As can be seen in Figure H] the distributions are non-normal in nature and often contain 
outliers away from the main peaks. Therefore, unlike the polarization angle distributions in 
Section 13. H none of these distributions are well characterized by a simple sample mean or 
a sample standard deviation. As robust descriptions of the distributions' central tendencies 
and width we use the samples' medians and median absolute deviations (MAD). The MAD 
of a set of measurements x is defined as the median value of the residuals, where the residuals 
are also calculated with respect to the sample median; that is 

MAD = median (|x -Xml) . (2) 

where Xm is the median value of the measurements x. For a normal distribution the MAD 
is significantly smaller than its standard deviation (a) with an expectation value of ct/1.48. 
However, given the small- number of comparison points in many of the clouds in Table [T] and 
that the P(850)/P(350)-distribution is not expected to be symmetric we report only the 
MADs there. 

The complete 3a data set contains 398 points with a median P(850)/P(350) value of 
1.9, MAD=0.7, and Xr=^- The x^-value implies that the distribution's width is intrinsic 
and not a result solely of data uncertainties. We reach the same conclusion examining the 
distributions for OMC-1, OMC-3, and DR21. (Medians are shown in Tabled = 12, 5, 
and 3.) 

We note that many of the peaks in the distributions of Figure H] are clearly different from 
the medians listed in Table [H this is mostly driven by some large outliers in the distribution. 
An alternate estimate of this peak is the value which minimizes the MAD. For the entire 
3a data-set in Figure IHa this alternative peak estimate is 1.5 with MAD=0.6. If the median 
in equation ([1]) is replaced with this peak then Xr = 6. These peaks, new MAD's, and Xr 
of OMC-1, OMC-3, and DR21 are given in Table [31 Using these peaks yields Xr > 1 ^'^^ 
three clouds and, therefore, does not change the conclusion that the distributions' widths 
are intrinsic and not a result solely of data uncertainties. Also, the difference between the 
median and the peak is less than the MADs in all cases (i.e., the whole data-set and the three 
specific clouds); none of the discussion in Section H] relies strongly on these precise values. 
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Fig. 4. — Histograms of the 850-to-350 /im polarization ratio for (a) all data in this work, (6) 
OMC-1, (c) OMC-3, and {d) DR21. Each plot contains different cuts on the data as shown 
in the legend of (c): all data (open histogram), only P > 3ap (light gray), P > Sap and 
\A(f)\ < 10° (dark gray). For reference vertical dotted lines are drawn at P(850)/P(350) = 1. 
For display purposes only, the insets in (a) and (6) magnify the lower-end of the distributions 
where the histograms peak. Note the larger plot in (d) contains all data in DR21 while the 
inset shows only data in DR21(Main). The bin-widths for (c) are 0.5, all others are 0.3. The 
total number of data points in each sample is given in Table [H 
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In using the polarization ratio to study grain alignment (e.g., Section 14. 2p we want 
to ensure that data at both wavelengths are sampling the same regions of the cloud, both 
along the line-of-sight (LOS) and across the plane-of-the-sky (POS). The latter criterion is 
simply met due to the fortuitous ability to beam-match the 350 and 850 /im data. If the 
emission sources are the same for radiation at both 350 and 850 /zm then the former criterion 
would also be met. Meeting this criterion is more difficult, but we try to limit its effect by 
choosing data with little-to-no position angle rotation between the two wavelengths. For this 
reason our analysis is often limited to data points where |A0| < 10°. Note that < 10° 
corresponds to data points with P > ScTp. 

To understand this particular data-cut, consider the case where the magnetic field 
changes its orientation along the line-of-sight (and within the cloud depths sampled by 
at least one wavelength). This may result in a wavelength-dependent change in both the 
polarization position angle (which follows the change of the field's projected orientation) 
and the polarization level (due to the field's changing inclination angle). In interpreting the 
polarization spectrum in terms of grain alignment physics, we wish to eliminate the changing 
inclination angle as the cause of any change in the polarization level (which can also result 
from other grain/cloud properties; see Section 1^^ . Since a changing field orientation must 
occur for any data with a wavelength-dependent angle, we can eliminate regions where this 
occurs by removing such data. While this does not ensure that points without wavelength- 
dependent angles arise from a single source it is unlikely, in the statistical sense, that the 
LOS field angles can change for many points in our large sample without an accompanying 
POS rotation. 

Table [1] shows the total number of data points satisfying both the P/o-p > 3 and |A0| < 
10° data-cuts in each cloud, along with their median polarization ratios and median absolute 
deviations. After making this second data cut the number of surviving data points drops 
to 141, only 35% of the 3a dataset. The resulting median ratio is P(850)/P(350) = 1.7 
with MAD=0.6 and Xr = H (see Figure la). The xl values for OMC-1, OMC-3, and DR21 
are also large using this data cut {Xr = 19, 5, and 3, respectively), again implying that the 
distributions' widths are not a result solely of data uncertainties. There are not as many 
outliers in the distributions after the |A0|-cut as their were in the 3(T-only cut. These do 
have some effect on the measured MADs (as discussed earlier in this section); Table E] reports 
the distribution peaks, MADs, and revised values of Xr for the |A0|-cut in the case where 
we have minimized the MADs. These small changes in the distribution widths still result in 
values of Xr ~ 1 meaning that our conclusion, that the widths are intrinsic and not a result 
solely of data uncertainties, still holds. 
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4. Discussion 



4.1. The Polarization Spectrum 



Fi^ureElshows an update d version of the polarization spectrum from I VaillancourtI (120021 ) 
and IVaillancourt et al.l ( 120081 ) . All the data in this figure satisfy the criteria P > 3ap and 
\A(j)\ < 10°. Here we plot the median value of DR21(Main), rather than DR21, as the data 
at other wavelengths (1300 /im) only cover that region of the cloud. While all the data in 
Figure E] show P(850) > P(350) we would draw the reader's attention to the range of ratios 
plotted in the distributions of Figure HI For example, the median ratio for all clouds in this 
work is plotted at P(850)/P(350) = 1.7 but has a relatively large MAD (0.6). 



The previous work comp aring Hertz and SCUBA-pol performed by IVaillancourtI (j2002[ ) 
and IVaillancourt et al.l (120081 ) were based on slightly different data sets than we use here. 
First, the ear lier results we r e obta ined usi ng the data prior to th e systematic re-analyse s 
performed by iDotson et al.l (120 lOt see also iKirby et al.l l2005l ) and iMatthews et al.l ( 120091 ). 
Second, those results attempted to match the Hertz and SCUBA beams by smoothing the 
SCUBA data to match Hertz's presumably larger beam size and re -sampled at a rate of 5 
arcseconds per pixel. Given the pointing accuracies of Hertz (4"-6"; iDowell et al.lll998[ ) and 
SCUBA (2")0 the 5" re-sampling was not unreasonable. However, as shown in Figure [21 the 
smoothing step was likely unnecessary. Despite these analysis differences the median results 
are in good agreement. For data satisfying P > 3(jp and |A0| < 10° the previous work 
found P(850)/P(350) medians an d standard deviatio ns of 1.4 ±0.6 and 1.7 ±2.7 for OMC- 
3 and DR 21 (Main), respectively ( IVaillancourtI |2002| ) . Here we find medians and MADs of 
1.6 ±0.6 and 1.6 ±0.4 for those two clouds. The large difference in the DR21(Main) standard 
deviations (=1.5 for the current work) likely lies in the fact t hat the earlier 8 5 0/350 data 
comparison used a different set of DR21 data. The a nalysis bv iMatthews et al.l (|2009[ ) used 
additional observations not available at the time of Vaillancourt ( j2002l ) and the analysis 
resulted in better rejection of noisy data. 



For W51, IVaillancourt et al.l ( l2008l ) found 1.8 ± 2.4 (median and standard deviation). 
This is consistent with the W51 results shown in Table [H 6.0 ±4.2 (median and MAD), but 
only because of the large deviations. These data are not plotted in Figure E] as only two data 
points survive the data cuts. 

The P(850) / P(350) ratio for OMC-1 is ca lculated in this work for t he first time. Figure[5] 
also includes P(450)/P(350) (solid triangle; IVaiUancourt et al.l l2008h and P(100)/P(350) 



^ http : //www . j ach . hawaii . edu/ JCMT/telescope/pointing/pointing_history . html 
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(open triangle; IVaillancourtI |2002| ) data points. The OMC-1 cloud is thus the only cloud 
in our data set which includes data both above and below the 350 /im minimum in the 
spectrum. 



4.2. Grain Alignment 



The original work by iHildebrand et al.l (Il999[ ) made it clear that the simplest model, 
isothermal dust populations all with the same polarization and/or alignment properties, can 
not reproduce a spectrum like that in Figure O In fact such a model yields a polarization 
spectrum independent of wavelength beyond 50 fim. In order to explain such a spectrum we 
consider that a number of physical mechanisms are responsible for the absolute polarization 
level observed in dust emission. Foremost among these are the efficiency with which dust 
grains become aligned with magnetic fields and variations in the inclination of that field 
to the line-of-sight. Ideally, our |A0| < 10° data cut (Section 13.21) has eliminated the field 
inclination as a variable in the observed spectrum, leaving alignment efficiency as the key 
variable. 

In order to generate wavelength-dependent polarization spectra the alignment efficien- 
cies must be correlated (or anti-correla ted) with changes in th e grains' emission. In order 
to explain spectra like those in Figure ISl lHildebrand et al.l (Il999[ ) considered simple emission 
laws of the form -F(i^) oc v^By{T)^ where v is the observed frequency, /3 the spectral index, 
and Bjj{T) is the Planck function at temperature T. For such models the required change 
in emission can take the form of differenc es in temperature, differences in the spectral in- 
dex, or a combination of the two (see also IVaillancourtll2002l 120071 ). Below we discuss two 
physical models of the ISM and molecular clouds, both of which result in grain popula- 
tions with the different alignment properties and temperatures/spectral-indices which lead 
to wavelength-dependent polarization spectra. 



Theoretical models of g rain alignment have a long history (see reviews by iLazarian 



2003 



2007 



Hildebrandl |1988[) with detailed observational tests possible only very recently 



2010 



Andersson et al. 



2011 



azarian. Goodrnan. &: Mverd 119971: iMatsumura fc Bastienll2009l : lAndersson fc Potter 



Matsumura et al. 



Hoang fc LazarianI l2008l . l2009a 



20111). In one of the most recent models , 



that of "radiative alignrn e nt torq ues" (RAT; ICho fc Lazarianll2005t iLazarian fc Hoang||2007 



b|) stellar and interstellar photons provide the necessary 



torque s to align the spin-axes of dust grains parallel to the local magnetic field. iBethell et al. 



( 120071 ) simulated a molecular cloud containing aspherical graphite and silicate gra ins with a 



typical interstellar grain-size distribution (i.e.. iMathis. Rumpl. &: Nordsiecki 119771 ) and radii 



of 0.005-0.5 /im. The RAT model results in alignment only of grain sizes larger than some 
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Fig. 5. — Polarization spectrum in several molecular clouds, normalized at 350 fim. The 
P(850)/P(350) data in this work are shown as dotted lines; note that OMC-1, OMC-3, and 
DR21(Main) have identical medians (see Tabled]). The solid circle represents the median of 
all data in this work. All data in this plot, including that shown here for the first time, meet 
the 3(Tp and A(f) criteria descr ibed in the te x t. Da ta at A < 350 fim and that for DR21(Main) 
at P(1300)/-P(350) are from IVaillancourtI ( 120021 ). All data used at 350 /im are from Hertz 
( iDotson et al.l |2010[ ) with the exceptior i of the OMC-1 point c omparing 350 and 450 fim 
(solid triangle) which is from SHARP (jVaillancourt et al.ll2008[ ). For clarity no error bars 
are shown here, but see Figure H] for the distributions. [^4 color version of this figure is 
available in the electronic version.] 
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cut-off, the exact value of wliicli is dependent on prop erties like tfie gas densi ty and radiation 
field and can vary throughout the simulated cloud (ICho fc Lazarianll2005l ). As the larger 
grains are more efficient emitters they reach cooler temperatures than the smaller grains 
in equilibrium. This yields an anti-correlation between grain temperature and alignment; 
the small warm grains are unaligned while the large cool grains are aligned. The resulting 
polarization spectrum rises from 100 to 400 /zm, with little variation at longer wavelengths. 
For the wavelengths of interest here they find P(850)/P(350) ~ 1.0-1.1. 



Draine fc Fraissd (120091 ) also present models composed of aspherical silicate grains and 



spherical graphite grains (here we discuss only their Model numbers 1 and 3). The grain-size 
distribution and the relative silicate-graphite mix are constrained by the observed interstellar 
extinction. Additionally, rather than model any physical alignment mechanism (such as 
RAT), their grain alignment is empirically constrained by the typical interstellar polarization 
spectrum spa nning near-opt i cal wa velengths (i.e., the "Serkowski law"). The result is similar 
to the work of iBethell et al.l (120071 ) in the sense that larger grains are both cooler and better 
aligned than the sm aller grains and produces a steep polarization spectrum in the 40 - 400 fim 
range. However, the lDraine fc Fraissd ( 120091 ) spectrum continues to rise beyond 400 /im such 
that P(850)/P(350) ~ 1.2-1.3. This long- wavelength behavior is a combination of a) the 
cooler silicate grains being aligned, whereas the warmer spherical graphite grains are not, 
and b) the shallower spectral indices of silicates compared to graphites. 

The median P(850)/P(350) values presented in this work (Fig. [S] and Table [T]) are 
clearly steeper than the model estimates just discussed (~ 1.7 for the "all clouds, P>3a, 
|A0|<1O°" sample). However, the data distributions are large and, therefore, cannot strongly 
rule-out either model. Additionally, the models are calculated using physic al conditions and 
constr aints which likely do not prevail in the real clouds studied here. The iDraine fc Fraisse 
(120091 ) model is constrained by data from the very low density ISM {Ay ^ a few) whereas 
our sample of clouds is flux-limited to s ome very bright, den se Galactic regions {Ay > 20). 
While modeling a "molecular cloud" , the iBethell et al.l (120071 ) model bathes the cloud rather 
uniformly in a typical interstellar radiation field which may be quite different from real clouds 
containing embedded stars. 



4.3. Embedded Sources 

The RAT mechanism predicts that grains exposed to stronger radiation sources will be 
more efficiently aligned. From this we might expect to see systematic trends in the polariza- 
tion with distance from stellar sources embedded in molecular clouds. Such a trend is hinted 



at in 60 /im polarization observations towards the W3A H II region (jSchleuning et al.ll2000l ). 
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The most prominent embedded sources in OMC-1 are a group of sources coincident with the 
BNKL intensity peak and the Trapezium stars in the visible Orion nebula (Fig. [6]). Using 
the MSX point-source cataloj^ ( Price et al. 2001 : Egan et al. 2003 ) we have also identified 
a number of embedded sources in the DR21 cloud (Fig. [7]). However, the proximity of the 
sources to each other, coupled with the fairly low spatial resolution of the polarization maps, 
does not allow a careful quantitative study of the strength of the polarization (at either 
wavelength) or the polarization ratio as a function of distance from these sources. 

A trend in polarization efficiency with distance from a radiation source also implies 
a correlation between the observed polarization and dust temperature. A careful measure 
of the dust temperature requires SED measurements over a wide range of wa. velengths, a 



task which is beyond the scope of the present work (e.g., I VaillancourtI |2002| ) . To some 
extent, one can consider the intensity or fiux density ratio, F(850)/F(350), as a proxy for 
the temperature. Figures [8]a and [8]& compare the intensity ratio in OMC-1 and DR21 to 
the polarization at both 350 and 850 ^m. The polarization in both clouds generally drops 
with increasing intensity ratio. If we interpret these ratios as color temperatures then the 
polarization increases with increasing temperature, as would be expected for grains aligned 
via RAT. However, we caution against over-interpretation of this result as defining a color- 
temperature is problematic in the dense clouds for at least two reasons. First, the range 
plotted in Figure |8] corresponds to unrealistically large temperatures; assuming (3 = 2 then 
T = 13K for the largest ratio F(850)/F(350) = 0.1 is reasonable but T > 100 K for ratios 
F(850)/F(350) < 0.033. Secondly, the ratio may also be the result of changes in grain 
emissivity (i.e., spectral index) and column density as well as the temperature. 

Using the intensity ratio as a proxy for temperature also has the advantage that it 
is independent of distance, allowing us to combine the relatively sparse data in individual 
clouds into a larger dataset. Combining the remaining data in clouds other than DR21 or 
OMC-1 results in Figure |H]c. The same trend of falling polarization is seen at both wave- 
lengths. We emphasize that we are comparing the polarization with the intensity ratio 
and are not discussing the "polarization-hole" effect which is often observ ed when compar- 
ing the polarization to absolute intensity at any given wavelength (e.g., 
Matthews. Fiege. fc Moriarty-Schievenll2002l ). 



Schleuning Il998 



ihttp : //irsa . ipac . caltech . edu/Missions/msx . html 
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Fig. 6. — Map of the polarization ratio, P(850)/P(350), in OMC-1. This map includes only 
P > Sdp data but includes points typically rejected by the |A0| criterion discussed in the 
text. Data with P(850)/P(350) > 3 are shown as saturated (black) pixels. Thin contours 
denote the 350 / xm intensity at levels of 1, 4, 8, 20, 40, and 80 % of the peak intensity (data 
from SHARC-2; IVaillancourt et al.ll2008[ l. The thick contour is drawn at P(850)/P(350) = 1. 
For reference we also show the four Trapezium stars of M42. [^4 color version of this figure 
is available in the electronic version.] 
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RA offset (arcmin) 

Fig. 7. — Same as Figure [6] but for DR21. Intensity data are from SHARC-2 (D. Dowell, 
private communication) Black stars indicate point sources identified by MSX (Section I4.3p . 
[^4 color version of this figure is available in the electronic version] 
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Fig. 8. — Polarization vs. intensity ratio in OMC-1 (a), DR21 (6), and all data except OMC- 
1 and DR21 (c). Here we plot only data meeting the P/(Jp and | A0| criteria described in the 
text. The lines are power-law fits to the data (solid line for 350 /xm, dashed line for 850 /xm); 
fits do not consider error bars. In (a) the fits exclude those points within 20" of OMC-l's 
central intensity peak (solid circles and solid triangles; see Section . 
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4.4. Intensity and Polarization Ratios 



One difficulty in using the absolute polarization values in Section 14.31 above is that 
the observed polarization magnitude is also a function of the parameters like the magnetic 
field's LOS inclination angle, grain cross-section, and turbulence, all of which may vary 
spatially across the cloud. The inclination angle effect is mitigated somewhat by our choice 
to limit the data set to those points with |A0| < 10° (see Section I3l2|) . These effects can 
be further mitigated by using the ratio P(850)/P(350). If the same grains are responsible 
for the polarized emission at both wavelen gths then those "polar ization reduction" factors 
effectively cancel in the polarization ratio ( iHildebrand et al.l 119991 ). 



Figures [6] and [7] show the spatial distributions of the polarization ratio in OMC-1 and 
DR21, respectively. Most of the mapped areas are characterized by polarizations which are 
larger at 850 /im than 350 /im. Notable exceptions are intensity peaks in both objects (BNKL 
at the origin of the OMC-1 map and DR21-0H(Main) at [+3'.'3, 0] in the DR21 map). The 
OM C-1 peak has also shown differences from the rest of the clou d in other polarization work 



l e.g.. lRao et al.lll998l : IVaillancourtll2002l : IVaillancourt et al.ll2008l ) so we will omit data within 
20" (one beam) of the peak in the analysis below. 

The most direct tests of the grain alignment models using submillimeter data require 
comparisons between the measured polarization ratio and the dust temperatures, spectral 
indices, and/or radiation environment of the aligned grains. A careful measure of those 
parameters requires SEDs measured o yer a wide range of wavelengths, a task which is beyond 



the scope of the present work (e.g., IVaillancourtI |2002| ). However, we can again use the 
intensity ratio, F(850)/F(350), as a proxy for the temperature or spectral index. Very 
different trends are observed when comparing this ratio to the polarization ratio in OMC-1 
and DR21 (Figures [9]a andP). The trend is an increase in P(850)/P(350) in OMC-1 but 
a decrease in DR21. As before the intensity ratio is independent of distance, allowing us 
to combine the remaining data in clouds other than DR21 or OMC-1 (Fig. [9]c). No strong 
trend between the intensity ratio and polarization ratio is observed. 

The large amount of scatter in these observations is not unexpected. The dense clouds 
studied here are certainly composed of multiple temperature components covering a wide 
range (e.g., ~ 20 - 80 K). Therefore, the intensity ratio is a measure not only of the 
components' dust temperatures, but also their different spectral indices and relative column 
densities. Presumably this effect is partly the cause of the large scatter observed in the data 
of Figure M 
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Fig. 9. — Polarization and intensity ratios in OMC-1 (a), DR21 (6), and all data except 
OMC-1 and DR21 (c). The upper x-axis scales cast the intensity-ratio as the standard 
spectral index defined through the relation F{i>) oc i/^"*"^. Here we plot only data meeting 
the P/(Tp and |A0| criteria described in the text. The lines are power-law fits to the data; in 
(a) the fits exclude those points within 20" of OMC-l's central intensity peak (solid circles; 
see Section 1131). 
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5. Summary 



We have comp iled all the spatially co incident data availabl e from the Hertz (350 /im; 
Dotson et al.ll2010[ ) and SCUBA (850 /im; iMatthews et aliboooh polarimeters. In order to 
facilitate this comparison the SCUBA-pol data have been spatially re-sampled to match the 
spatial locations of the Hertz data; the two instruments have comparable spatial resolution 
(20"). We find a total of 1699 individual locations which can be compared within 17 different 
Galactic clouds; all the polarization and intensity data at 350 and 850 fim are given in 
electronic format in Table 2. Of these data 1124 points yield non-zero polarizations at 
both wavelengths when corrected for the positive polarization bias. Additionally, of the 
bias-corrected data, 398 points satisfy the criterion P/cXp > 3 at both wavelengths, and 141 
points satisfy the additional criterion |0(85O) — 0(350) | < 10°. Complete polarization and 
intensity maps for all clouds are given in the Appendix. 

We have investigated the change in polarization angle from 850-to-350 /xm. The angle 
differences exhibit a wide distribution indicating that, in some regions, there is a real angle 
rotation within the measurement uncertainties. However, the angle distributions are centered 
about 0(350) ~ 0(850). These conclusions hold globally for the entire -P/cXp > 3 data set 
in this work and individually for the P/(Jp > 3 data sets in the clouds OMC-1, OMC-3, 
and DR21. Due to the limited number of data points in other individual clouds we have 
made no attempt to study the angle distributions nor the point-by-point angle agreement 
between wavelengths; therefore these conclusions do not necessarily extend to the other 
clouds included herein. 

We have also examined the 850-to-350 /im polarization ratio, P(850)/P(350), on a point- 
by-point basis at every spatial location in each of our sampled clouds. From this work we 
establish a genuine trend towards higher polarization at 850 fim than 350 fim with a median 
polarization ratio of P(850)/P(350) = 1.7 (and a median absolute deviation of 0.6). These 
values cover all clouds in our dataset and points which satisfy the criteria P/ap> 3 and also 
85 0) — 0(350)1 < 10°. This trend is consistent w ith previous work (e.g.. iHildebrand et al. 



19991 : IVaillancourtI l2002l : IVaillancourt et al.l l2008l ) and is best explained by models which 



require mixtures of dust grains with different physical properties (i.e., temperatures and 
spectral indices) an d dif ferent alignment e fficien cies. The dust and alignment models of 
Bethell et al.l (120071 ) and iDraine &: Fraissd (120091 ) predict increasing values for the polariza- 
tion from 350-to-850 fim (~ 1.0-1.3) but not of the same magnitude as observed here. This 
is most likely due to the fact that their models use very different physical conditions than 
prevail in our sample of bright, dense molecular clouds. 

We find a trend in which the 350 and 850 /im polarizations tend to fall as the 850-to-350 
fim intensity ratio increases. If we interpret this ratio as a color temperature then these data 
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are consistent with a key prediction of radiative alignment torques. That is, grains which 
arc more exposed to radiation sources, and are thus warmer, are more efficiently aligned. 
However, we caution that the assignment of temperature to a two-wavelength intensity ratio 
is not robust. No clear trends arc observed when the polarization ratios are compared to the 
intensity ratios on a point-by-point basis. Better tests require work at additional wavelengths 
in order to produce SEDS from which more accurate dust temperatures can be extracted. 

This work would not have been possible without the dedicated support of the staff at 
the Caltech Submillimeter Observatory and James Clerk Maxwell Telescope and the ex- 
tended team of dedicated scientists who toiled to collect and analyze data for both Hertz 
and SCUBA-pol. We would hke to thank Roger Hildebrand, Giles Novak, and B-G An- 
dersson for comments on an early draft of this paper. This research has made use of the 
NASA/IPAC Infrared Science Archive, which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with the National Aeronautics and Space 
Administration . 
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Table 3. Polarization Ratio Distributions 



Source 


data 


satisfying 


P > 3o-p 


...also 


satisfying 


A(j)\ < 10° 


Peak 


MAD 




Peak 


MAD 




OMC-1 


1.3 


0.5 


13.4 


1.4 


0.5 


20.8 


OMC-3 


1.4 


0.5 


2.9 


1.2 


0.3 


3.7 


DR21 


1.4 


0.4 


2.5 


1.8 


0.3 


5.2 


DR21(Main) 


1.3 


0.3 


2.3 


1.8 


0.3 


6.6 


All 


1.5 


0.6 


6.0 


1.4 


0.5 


10.9 



Note. — The peak value and the median absolute deviation (MAD) of the 
polarization ratio (P[850]/-P[350]) distributions which minimize the MAD 
(eq. [2]). Also shown are Xr values as calculated from equation ([1]); see 
Section [321 The columns labeled "P > 3cTp" and "also \A4>\ < 10°" are 
defined as in Tabled) 
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A. Polarization Maps 



Figures [T0] - [26] present grayscale/contour maps of intensity along with polarization vec- 
tors at 350 /im (blue lines) and 850 /im (red lines). Only vectors with P > ScXp and cxp < 4% 
are drawn on the maps; the latter criterion is an aesthetic choice to remove points with atyp- 
ically large polarizations. Vectors spaced more closely than the nominal 17'.'8 Hertz pixel 
pitch are a result of moving the Hertz array in sub-pixel steps. Additionally, some SCUBA- 
pol vectors in regions without Hertz coverage have not been plotted. Each map includes a 
scale-bar for determining the absolute polarization levels; the scale-bars differ for each map 
but are the s ame length for each w avelength within a map. Map center coordinates are given 
in Table 1 of IPotson et~aD toid ). 



As indicated in the figure captions most intensity maps are from SCUBA at 850 /im 
(jPi Francesco et al.ll2008l ). For aesthetic reasons some maps use 350 /xm intensity data from 
Hertz, wh ile OMC-1 and DR 21 use 350 yum intensity data from the SHARC-2 camera at 
the CSO (jPowell et al.ll2003l ). The lower-right corner of each map includes a gray-circle 
indicating the 20" effective beam-size of Hertz and SCUBA-pol. Note that this beam is for 
the polarization data only, not the intensity maps; SCUBA intensity maps typically have 
19" resolution. Hertz intensity maps typically have 28" resolution, and SHARC-2 intensity 
maps typically have 10" resolution. 
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Fig. 10. — W3. The intensity map is from Hertz at 350 /xm with contours drawn at 10, 20, 
30, 90% of the peak intensity. Other key map features for this and subsequent figures are 
described in the text. 
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Fig. 11. — NGC 1333. The intensity map is from SCUBA with contours drawn at 2, 5, 10, 
20, 40, 60, and 80% of the peak 850 fim intensity (which occurs in IRAS4A). 
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Fig. 12.— OMC-1. The intensity map is from SHARC-2 f IVaillancourt et all l2008l ) with 
contours drawn at 1, 2, 4, 6, 8, 10, 20, 40, 60, and 80 % of the peak 350 /xm intensity. The 
gray circle in the lower-right indicates the 20" effective beam-size of Hertz and SCUBA-pol, 
not the SHARC-2 intensity data which has a beam size of ~10". 
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Fig. 13. — OMC-2. The intensity map is from SCUBA with contours drawn at 5, 10, 20, 30, 
90% of the peak 850 /xm intensity. 
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Fig. 14.— OMC-3. The SCUBA intensity is shown with contours at 10, 20, 30, 40, 60, and 
80% of the peak 850 /im intensity. (Note that some SCUBA-pol vectors in the center of the 
map have been removed where there was no Hertz-coverage of this cloud.) 
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Fig. 15. — NGC2024. The intensity map is from SCUBA with contours drawn at 5, 10, 20, 
30, 90% of the peak 850 /im intensity. 
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Fig. 16. — NGC2068 LBS 10. The intensity map is from SCUBA with cont ours drawn at 
10, 20 , 30, 90% of the peak 850 /im intensity. Cloud core labels are from [Mitchell et al. 
(|200l[ ). 
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Fig. 17. — NGC2068 LBS 17. The intensity map is from SCUBA with cont ours drawn at 
10, 20 , 30, 90% of the peak 850 /im intensity. Cloud core labels are from [Mitchell et al. 
(|200l[ ). 
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Fig. 18. — NGC2071. The intensity map is from Hertz with contours drawn at 20, 30, 
90% of the peak 350 /im intensity. 
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Fig. 19. — MonR2. The intensity map is from SCUBA with contours drawn at 10, 20, 30, 
90% of the peak 850 fim intensity. 
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Fig. 20.— MonOBl IRAS 12 (IRAS 06382+0939). The intensity map is from SCUBA with 
cont ours drawn at 10, 20, 30 , 90% of the peak 850 /xm intensity. Cloud core labels are 
from IWoIf-Chase et al.l J2003h 
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Fig. 21. — pOph A. The intensity map is from SCUBA with contours drawn at 10, 20, 30, 
90% of the peak 850 pm intensity. 
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Fig. 22. — IRAS 16293—2422. The intensity map is from Hertz with contours drawn at 10, 
20, 30, 90% of the peak 350 /im intensity. 
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Fig. 23. — NGC6334A. The intensity map is from SCUBA with contours drawn at 5, 10, 
20, 30, 90% of the peak 850 /im intensity. 
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Fig. 24.— W49A. The SCUBA intensity is shown with contours at 10, 20, 30, 90% of 
the peak 850 fim intensity. 
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Fig. 25.— W51 A (G49.5-0.4). The Hertz intensity is shown with contours at 10, 20, 30, 
90% of the peak 350 /im intensity. 
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Fig. 26. — DR21. The 350 fim intensity is shown in grayscale with contours drawn at 1, 2, 
5, 10, 20, 40, and 80 % of the peak. Intensity data are from the SHARC-2 camera with a 
spatial resolution of ~ 10" (C. D. Dowell, private communication). The gray circle in the 
lower-right indicates the 20" effective beam-size of Hertz and SCUBA-pol, not the SHARC-2 
intensity data. The intensity peak is DR21 OH (Main) and the southern peak at the map 
origin is DR21 (Main). 
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